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(Abstract) An Uzawa-SOR method is presented for solving fuzzy linear system whose coefficient matrix is crisp and the right-hand 
side column is arbitrary fuzzy number vector. The convergence is analyzed and numerical example is given to illustrate the procedure. 
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1. INTRODUCTION 

Fuzzy linear system (FLS) has many applications in control 
problems, information, statistics, engineering, economics, finance 
and even social sciences. In the 1990s, Buckley et al. [8-10] 
investigated them in series. Rao and Chen [16] consider the 
numerical solution of FLSs in engineering analysis. Friedman et 
al. [12] consider a FLS as follows, 
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Definition 2.1. A fuzzy number vector X = (x l ,X 2 ," 
given by 

x t =(x i (r),x i (r)), \<i<n,0<r<\, 
is called a solution of FLS (1.1) if 
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By (2.1), Friedman et al. [12] extend FLS (1.1) to a In x In 
crisp linear system 

SX = Y (2.2) 



where S = (s kl ), S kl are determined as follows 



Sij =a ip 



S i+n,j+n a ij % - 



^i,j+n i+nj 



1 < j, j < n, 



fl ,A+V2+-+v„ = y^ 

where the coefficient matrix A = (a^ ) is a crisp matrix and 
y = (j.) is a fuzzy vector, 1 < i, j < n. They suggest a 
model to solve this kind of fuzzy system. Based on it, numbers of 
numerical methods [1-7,11,13-15,17-19] have been presented 
for FLS (1.1). In this paper, we provide an iterative method 
named Uzawa-SOR (cf. [20] ) for solving FLS (1.1). 

The paper is organized as follows. In Section 2,we give some 
preliminaries for FLS (1.1). In Section 3, we propose the 
Uzawa-SOR method with the convergence theorem. An 
illustrative example is provided in Section 4 and the conclusion is 
drawn in Section 5. 

2. PRELIMINARIES 

Following [12], a fuzzy number is defined as (u(r),u(r)) , 
< r < 1 , which satisfies, 

• u(r) is a bounded left continuous nondecreasing function 
over [0,1], 

• u(r) is a bounded left continuous nonincreasing function 
over [0, 1], 

• u(r) < u(r)), < r < 1 . 

To define a solution to (1.1) we should recall the arithmetic 
operations of arbitrary fuzzy numbers X = (x(r), x(r)) , 
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(2.1) 



a, y >0 
a tj <0 



and any S kl which is not determined by the above items is 
zero, 1 < k,l < 2fl, and 
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From [12], we know that S has the following structure 

S { S 2 
_S 2 5j_ 

where 5 p 5 2 >0, A = S { — S 2 , and (2.2) can be 
rewritten as follows 

\S,X-S 2 X=Y, 

- - (2.3) 
[S l X-S 2 X=Y, 

where 
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The following theorem indicates when FLS (1.1) has a 
unique solution. 

Theorem 2.2 [12]. The matrix S is nonsingular if and only if the 
matrices A = Si - S 2 and Si + S 2 are both nonsingular. 

Under the conditions of Theorem 2.2, the solution X = S' l Y 
of ( 1 . 1) is thus unique but may still not be an appropriate fuzzy 
vector. By Theorem 2 of [12], we know that S~ l has the same 
structure like S, i.e. 

, fa T 2 

T T 

_ 2 1 \_ 

The following result provides a sufficient condition for the 
unique solution to be a fuzzy vector. 

Theorem 2.3 [12]. The unique solution X of (2.2) is a fuzzy 
vector for arbitrary fuzzy vector Y, ifS' 1 is nonnegative. 

Restricting the discussion to triangular fuzzy numbers, i.e. 

y ,{r), y, (r) and consequently x ( (r), 3c ; (r) are all linear 

functions of r, and having calculated X which solves (2.2), we 
can define the fuzzy solution to the original system given by 
(1.1) as follows. 

Definition 2.4. LetZ = { (x, (V),-^- (r)), 1 < i < n} denote 
the unique solution of (2.2). The fuzzy number vector 
U = {(w ; (r), w ; (r)), 1 < i < n) defined by 

(r) = min{ x. (r), x.(r), x t (1), x t (1) } , 
u f (r) = maxiXfirXXiirXXfO-XXfO)} 
is called the fuzzy solution of SX = Y. If (x t (r), x ( (r)), 
\<i<n are all fuzzy numbers then u t (r) = x t (r), 
M ( (r) = 3c,. (r), 1 < i < n and U is called a strong fuzzy 
solution, otherwise, U is called a weak fuzzy solution. 
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3. THE UZAWA-SOR METHOD 

For the case S is nonsingular, without loss of generality, assume 
that s u > 0, i = 1,2, ■ • • , n, and S 1 = D l - L x -U l , where 

= diag(^ |( ), i = 1,2, L x and £/ t are strictly lower 

and upper triangular matrices. Let < CO < 2 be a relaxation 
parameter. Then for the first equation of (2.3), we take the 
following SOR iterative form: 

(D, - coL^X^ = [(1 - a)D l + aUJX, + co(S 2 X k + Y), (3.1) 

and for the second equation of (2.3), we can take the iterative 
scheme as follows, 

X k+I =X k + T(S 2 X M - S,X k - Y), (3.2) 

where T is a real parameter. 

Then we get the Uzawa-SOR method for FLS (1.1) as the 
following algorithm. 

Algorithm 3.1. (Uzawa-SOR method) Given initial vectors 
X_ ,X Q £ R", a relaxation factor < CO < 2 and a real 

parameter T. For k = 0, 1, 2, • • • , the following iterative 
scheme is taken, 

X k+l = (D, - coh x y x [(1 - a>)A + coU, ]X k 

+ co(D l -coL i y l (S 2 X k +Y_), (3.3) 
%k+i = Xk + T (SiX-k+i ~^iX k -Y). 
We have the following convergence theorem. 
Theorem 3.2. If S 1 is symmetric positive definite, let X m 
and X M denote the smallest and the largest eigenvalues of S x , 
then for < CO < 2 and < T < ~r~ , the Uzawa-SOR 
method (3.3) is convergent. The optimal parameter T is 

V = arg min 0< r< ^_ { 1 1 - tA m 1, 1 1 - rA m I } . 

Proof. Because Si is symmetric positive definite, for 
< CO < 2 , the SOR iteration (3.1) is convergent. For 
iteration (3.2), we can rewrite it as 

^=(7-^ + 1(5^-7). 

The iteration matrix of (3.2) is I — tS 1 . Suppose A be an 
arbitrary eigenvalue of I — rS l and z be the corresponding 
eigenvector, we have 

(I-tS 1 )z = Az, 

that is 



Let X s be an arbitrary eigenvalue of Su thus A = 1 — zA Si . 
Therefore 
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I X I «■ \1-tX s kl 



<=> 0< r < . 



\1-tX s < 1 
1 1 - tX s > -1 



4, 

Then we see that the iterative scheme (3.2) is convergent if 
< r < -j^ , and the optimal T is arg min o< t< ^_ { I X I } , i.e. 



argmin 0<r< ^{ll-r/l M l,ll-r/l m l}. 



□ 



4. NUMERICAL EXAMPLE 
Example. Consider 2X2 fuzzy linear system 

fx, -x 2 =(r,2- r), 

[x, + 3x 2 =(4 + r,7-2r). 

The exact solution is 

fx, = (x, (r), 3Cj (r)) = (1 .375 + 0.625 r, 2.875 - 0.875 r), 
|x 2 = (x 2 (r), 3c 2 (r)) = (0.875 + 0.125 r, 1.375 - 0.375 r). 

Direct calculations yield X m =0.8453, X M =3.1547, 

r opt =0.5. 

With C0=\,T = r opt = 0.5 and X Q = [0,0,0,0] T , after 
17 iterations, the numerical solution is 

"1.3748 + 0.625 lr" 
0.8751 + 0.1250r 
2.8748 -0.8749r 
u 1.3753-0.3751r_ 

i.e. 

fx, = (1.3748 + 0.625 lr, 2.8748 - 0.8749r), 
[x 2 = (0.8751 + 0.1250r, 1.3753 - 0.375 lr). 

With co = 0.95, r = r opt = 0.5 and X Q = [0,0,0,0]' 

after 17 iterations, the numerical solution is 

"1.3750 + 0.6250r" 
0.8750 + 0.1250r 
2.8748 -0.8749r 
1.3753 -0.3751r 



X 



1,0.5,17 



X 



0.95,0.5,17 



i.e. 



x, = (1 .3750 + 0.6250r, 2.8748 - 0.8749r), 
x 2 = (0.8750 + 0.1250r, 1.3753 - 0.375 lr). 



5. CONCLUSION 

We present an Uzawa-SOR iterative method for nx n fuzzy 
linear system. If the proposed matrix S by Friedman et al. [12] 
is nonsingular, then for any initial vector X Q , the Uzawa-SOR 
iteration will converge to the unique solution of SX = Y. The 
numerical example shows that the method is effective, and with 
appropriate parameters CO and T, it is expected to achieve 
higher accuracy or faster convergence rate than the known 
methods which have one or no parameter. As analyzed in 
Section 3, it is just obtained the convergence range for the 
symmetric positive definite case and the optimal parameter X . 
The convergence for general case and the optimal CO need to 
be further investigated. 
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By taking various CO and r, we can get higher accuracy or 
faster convergence rate. 

IJICS Volume 1 , Issue 2, May 201 2 PP. 30-33 www.iji-cs.org © Science and Engineering Publishing Company 

-32- 



International Journal of Information and Computer Science 

Syst. 11 (2009) 308-313. 
[14] H.-K. Liu, T.-T. Chang and B. Wu, Numerical methods for 

solving fuzzy system of linear equations, ICIC Express Letters 4 

(2010) 25-30. 

[15] S.-X. Miao, B. Zheng and K. Wang, Block SOR methods for 
fuzzy linear systems, J. Appl. Math. Comput. 26 (2008) 201-218. 

[16] S.S. Rao and L. Chen, Numerical solution of fuzzy linear 
equations in engineering analysis, Internat. J. Numer. Methods 
Engrg. 42 (1998) 829-846. 

[17] K. Wang and B. Zheng, Symmetric successive overrelaxation 



IJICS 

methods for fuzzy linear systems, Appl. Math. Comput. 175 
(2006) 891-901. 

[18] K. Wang and B. Zheng, Block iterative methods for fuzzy linear 
systems, J. Appl. Math. Comput. 25 (2007) 119-136. 

[19] J.-F. Yin and K. Wang, Splitting iterative methods for fuzzy 
system of linear equations, Comput. Math. Model. 20 (2009) 
326-335. 

[20] J. Zhang and J. Shang, A class of Uzawa-SOR methods for saddle 
point problems, Appl. Math. Comput. 216 (2010) 2163-2168. 



IJICS Volume 1 , Issue 2, May 201 2 PP. 30-33 www.iji-cs.org © Science and Engineering Publishing Company 

-33- 



